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I. 


INTRODUCTION 


In  recent  years  the  geometrical  theory  of  diffraction  [1]  has 
become  one  of  the  more  useful  methods  for  finding  solutions  to  antenna 
and  scattering  problems.  The  geometrical  theory  of  diffraction  is, 
however,  limited  to  geometries  for  which  the  diffraction  coefficient 
is  known. 

Quite  recently  a technique  has  been  developed  for  extending  the 
GTD  by  the  method  of  moments  [2].  In  this  technique  the  diffraction 
coefficient  is  treated  as  an  unknown,  thereby  permitting  a larger  class 
of  problems  to  be  treated  with  the  GTD.  In  the  case  of  planar  geo- 
metries such  as  the  wedge,  this  procedure  has  been  improved  by  Sahalos 
and  Thiele  [3]  using  a three  term  representation  of  the  diffraction 
coefficient.  Using  this  improved  technique,  it  has  been  shown  that 
one  can  find  the  current  near  a discontinuity  in  a surface,  such  as 
near  the  edge  of  a wedge,  for  an  arbitrary  incidence  angle. 

Previously,  the  extended  GTD  technique  has  been  applied  to  2- 
dimensional  problems.  In  this  report  we  will  use  the  method  to  treat 
a 3-dimensional  geometry,  the  three  sided  pyramid.  The  reason  for 
treating  this  particular  geometry  is  to  demonstrate  the  applicability 
of  the  extended  GTD  technique  to  3-dimensional  geometries.  As  such, 
it  represents  a first  step  in  applying  the  GTD-moment  method  technique 
to  3-dimensional  problems.  Other  geometries  and  modifications  of  the 
method  to  be  described  in  this  report  will  be  considered  subsequently. 

For  the  problem  of  scattering  by  a pyramid  as  treated  in  this 
report,  each  face  of  the  pyramid  is  composed  of  two  regions,  a GTD 
region  and  a moment  method  region  near  the  edges.  The  current  distri- 
bution in  the  GTD  region  is  found  by  solving  the  3-dimensional  wedge 
diffraction  problem  once  for  each  wedge,  or  a total  of  six  times. 

Next,  the  current  near  the  edge  of  each  wedge  is  found  using  the  moment 
method.  However,  the  currents  obtained  on  the  faces  of  the  pyramid 
do  not  include  that  current  caused  by  tip  diffraction.  The  effect  of 
this  current  on  the  far  zone  scattered  field  may  be  added  post  facto 
using  a formula  by  Keller,  et  al  [4].  While  the  effect  of  tip  dif- 
fraction is  not  large  [5],  its  inclusion  does  bring  the  theoretical 
results  into  closer  agreement  with  measurements. 

In  the  work  done  here,  only  first  order  diffraction  is  considered 
and  this  is  found  to  be  adequate  for  determining  the  far  zone  field. 

On  the  other  hand,  if  one  was  interested  in  the  near  zone  field,  it 
would  be  imperative  that  at  least  second  order  diffraction  be  con- 
sidered. 

Section  III  presents  the  extended  GTD  formulation  for  the  pyra- 
mid while  Section  IV  gives  results  for  the  bistatic  radar  cross  section 
as  well  as  the  currents  on  the  cone.  Some  concluding  remarks  are  made 
in  Section  V. 
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II.  FORMULATION 


Consider  the  geometry  of  the  three  dimensional  wedge  diffraction 
problem  (Figure  1).  The  field  at  any  point  in  space  will  be  an  expres- 
sion of  the  form 


U(r,(ji)  = V(r,t+4>')  + V(r, ) (1) 


where  V(r,(ji+4>')  defines  the  field  due  to  the  incident  wave  and  V(r,4>-'}'') 
the  field  due  to  the  reflected  wave. 

If  we  want  to  have  an  expression  for  the  geometrical  optics  field 
and  the  diffracted  field  separately,  then  by  the  help  of  the  uniform 
theory  of  the  diffraction  [1]  we  can  see  that: 

(2) 

where  in  the  general  situation  the  function  will  be 
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Vg(L,(|)*)  = I_„(L,(j>^)  + I+.„(L,i|)^)  . (3) 
where  (j."*"  = and  (>>"  = 

The  parameter  L is  not  the  distance  of  the  observation  point  but  a 
function  of  this  distance  depending  upon  the  nature  of  the  incident 
wave. 


Now  U-CL.iC-)  has  been  given  by  Kouyoumjian  as  follows: 


j(0L+Tr/4) 


(8La)l/2 


where  B»2it/x  and  a*l+cos(i^~-2nirN). 

N is  an  integer,  positive,  negative  or  zero,  whichever  best  satisfies 
the  equations 

2nitN-^^  ■ -IT  , for  (5) 


1 

1 

I 

I 

1 

1 

.1 


and 


-■} 


2 


DIFFRACTION  POINT 


Figure  1.  Geometry  for  three-dimensional  wedge 
diffraction  problem. 
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(6) 


2miN-(|)'  =+Ti,  for  1+^ 


Analyzing  the  incident  and  the  diffracted  field  in  terms  of  the  polar 
components  at  the  observation  point  we  will  have 
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where 


A(s)e-jes  (7) 

smBp 


= Vgd.D  * Vgd,/) 


(8) 


and  the  spatial  attenuation  factor  A(s)  is  defined  as 


for  plane,  cylindrical  and  conical  wave 
incidence 

for  spherical  wave  incidence 


(9) 


and  the  distance  parameter  L 
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L = 
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s sin^Bo 

r'  r 
r+r* 

S'  s sin^Bp 


is  defined  as 

plane  wave 
cylindrical  wave 

conical  and  spherical  wave 


(10) 


Let  us  now  consider  a triangular  pyramid  (Figure  2).  If  we  want 
to  study  the  diffraction  we  can  easily  see  that  the  diffracted  field 
at  any  point  will  be  the  summation  of  six  diffracted  fields  given  by 
Equation  (7).  Usually  we  need  to  find  the  radar  cross  section,  or  the 
radiation  of  antennas  which  are  on  the  body.  So,  it  is  helpful  to  con- 
sider the  current  on  the  body.  This  current  gives  the  necessary  informa 
tion  for  the  electromagnetic  properties  of  the  body. 

Because  of  the  incident,  reflected  and  diffracted  fields,  we 
suppose  that  on  the  body  we  have  surface  currents  normal  to  the  surface 
fields.  The  fields  on  the  surface  of  the  body  can  be  computed  by  the 
help  of  the  aoove  equations.  By  knowing  the  surface  current  we  can  find 


4 


the  field  at  any  point  in  space.  The  current  on  the  surface  is  depend- 
ent only  on  the  field  and  will  be  the  vector  summation  of  currents 
because  of  incident,  reflected  and  diffracted  rays.  In  the  problem  of 
a pyramid,  for  example,  the  current  on  the  point  M (Figure  3)  will  be 
the  suinmation  of  four  currents.  ll,  are  currents  due  to  the 

diffraction  on  the  edges  1,  2 and  4 correspondingly  and  + is  the 
current  due  to  the  incident  and  reflected  fields. 

All  the  above  expressions  can  give  us  is  an  approximation  for  the 
electromagnetic  prope-^ties  of  the  pyramid  and  the  results  could  be 
sufficient  if  the  uniform  theory  of  the  diffraction  was  valid  for  points 
near  the  edges.  Since  it  is  not,  a MM-GTD  formulation  is  necessary. 


III.  FORMULATION  OF  MM-GTD  SOLUTION 

A.  Three  Dimensional  Wedge  Diffraction 

Let  us  consider  again  (Figure  4)  the  geometry  for  the  three  dimen- 
sional wedge  diffraction  problem.  To  find  the  surface  current  on  the 
walls  we  can  use  the  magnetic  field  integral  equations.  By  the  assum- 
tion  that  the  incident  wave  is  a plane  wave,  the  current  parallel  and 
normal  to  the  edge  has  been  found  by  Sahalos  and  Thiele  [2].  The  problem 
is  reduced  to  the  two  dimensional  one  and  the  method  of  moments  along 
with  the  GTD  as  used  by  Burnside,  et  al  [3]  gives  very  good  results. 

The  current  on  the  XiY-j  plane  is  expressed  by  orthogonal  pulses 
in  the  MM  region  and  by  GTD  elsewhere.  So  the  jj(|  and  Jyj  current 
parallel  to  (Xi.Y-j)  will  be  given  by  the  solution  of  linear  equations 
as  below [2]: 


where  a^,  is  the  weight  of  the  n^h  pulses  in  MM  region  and  D*^  is  the  dif 
fraction  coefficient  in  the  GTD  region. 
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The  current  is: 


n=l 


1 „-jkcos6«. 

jUjr  , I Dj  e 

i=-l  /F 


MM- region 


GTD-region 


where  t is  either  X-]  or  Yi. 

The  and  currents  can  be  expressed  [2]  as  a combination  of  a^, 
Df,  and  they  have  been  given  also  by  Sahalos  and  Thiele.  The  parameters 
I ' , kp,  IpiQ.  Iq  and  kp  are  dependent  on  the  geometry  of  the  v;edge 
and  their  expressions  are  in  Reference  [3]. 


Figure  4.  Geometry  showing  6TD  and  MM  regions 
on  both  sides  of  a wedge. 


Figure  5.  Geometry  showing  two  faces  (ABC  and  ACD) 
of  the  pyramid  for  far  zone  scattering. 

Suppose  that  we  want  to  find  the  far  zone  scattered  field  in  the 
direction  (41  ,e^)  due  only  to  the  surface  ABCD  which  is  composed  of  two 
triangular  plates  as  shown  in  Figure  5.  We  will  use  the  field  expression 


r 

>^0 


(ABCD) 


(r.  X Jl)e 


jkr-r. 


where  r^  is  the  distance  between  the  origin  point  A and  the  observa- 
tion point,  Jp  is  the  current  on  the  wall  and  k is  the  wavelength  number. 
The  current  Jj*  on  the  two  faces  ABC  and  ADC  can  be  exp<'essed  as: 

g-jkzcosB 

'Jp  * [J2r(x')z  Jy, (x)cosY  X - Jy. (x')slnY  y]  e^JkzcosB  ^ 
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The  vector  Tj  can  be  expressed  in  the  three  components  parallel  to  x, 
y,z  by  the  help  of  So, 


r^  = sine^cosiii^  x + sine^siniji^  y + cose^  z (15) 


] 


= - siny 

(19) 

A^  ( 415 » Oj ) 

- + A2((t>jiY,0j) 

(20) 

can  give  the  expression  of 

The  vector  expressed  in  another  cartesian  system  A (xQyQZg) 

will  be  given  by  using  the  Euler's  angles.  So 


(21) 


where  Rq  is  the  distance  AAq  and  K is  the  matrix  of  the  directional 
cosines  of  the  axis  x,y,z  to  the  axis  Xo,yo»Zo-  field  will  be: 


/hJ\ 

*0 

0 -cose^  sine^sin((i^ 

K\ 

j(k-R  -kr  ) 

J'o 

g 

2roX 

-cose^  0 -sine^cos(fi^ 

yj 

-sin05sin((i^  sine^cos^^  0 

(22) 

At  this  point  we  have  only  the  scattered  field  from  two  triangular 
plates  by  using  the  diffraction  from  the  common  edge.  In  the  next 
section  we  will  use  the  above  expressions  to  obtain  the  diffracted  field 
from  a pyramid. 


B.  Diffraction  By  a Pyramid 

A pyramid  is  a combination  of  six  edges  and  since  the  length  of 
the  edges  are  more  than  the  wavelength  the  problem  of  the  diffraction 
can  be  separated  into  six  problems  of  wedge  diffraction.  For  the  six 
edges  (see  Figure  6)  the  scattered  field  will  be  the  summation  of  six 
expressions  like  Equation  (22). 
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Figure  6.  Six  diffracting  edges  of  the  triangular 
pyramid. 


First,  it  is  important  to  define  the  shadow  region  of  the  pyramid. 
With  the  help  of  Figure  7 we  can  see  for  example  that 

if  41 1 2 ^ ^ the  face  12  is  shadowed 

if  < i»-A^  the  face  13  is  shadowed 

A more  general  definition  of  the  shadow  region  can  be  obtained  in  the 
following  manner.  We  take  the  horizontal  plane  (456)  and  starting  from 
the  direction  of  the  incident  wave  we  measure  the  angles  4)1,412  from  the 
edges  1 and  2.  We  define  the  maximum  of  the  Ratios: 


12 


tane^'^  tanei^ 

1 ’ ” 2 
tane'  tane^ 


tane^^ 

tane^ 


\ 


If  any  of  the  is  not  in  the  distance  (-1,1)  then  we  are  not  in  the 
shadow  region  except  the  face  (456).  If  all  t,  are  in  the  distance 
(-1,1)  then  w define 


i 


1 


cos<>5  = ax(T^) 


and  we  construct  the  following  table: 


Illuminated  1 

Angle 

edges 

faces 

1 

0<ij>g£4)i 

1+3+6 

(13)  for  e^"<iT/2 

1+3+4+5+6 

(13)+(456)  for  0^">ir/2 

2 

^2 

1+2+3 

(12)+(13)  for  e’'’<iT/2 

1+2+3+4+5+6 

(12)+(13)+(456)  for  0’'’>7r/2 

3 

1+2+3 

(12)+(13)+(23)  for 

From  this  table  we  can  determine  if  the  wave  is  striking  face  (23)  or 
face  (12). 

The  above  table  is  useful  for  computations  because  with  the 
described  conditions  we  can  easily  program  the  diffracted  and  G.O. 
fields.  The  incident  ^ield  Hi  (Figure  4)  is  normal  to  the  vector  R 
and  the  definition  of  e"  gives  the  value  of  with  the  help  of  the 
formula: 


l^i.^Hj  , cos"^(cot0^cot8^) 


(23) 


In  our  problem  of  finding  the  radar  cross  section  the  diffraction 
from  the  tips  is  negligible  in  comparison  to  the  total  scattered  field. 
So,  except  for  the  case  of  axial  incidence,  by  the  above  work  we  can 
have  a good  approximation  for  the  radar  cross  section  without  considering 
tip  diffraction.  In  a future  report  we  will  show  a method  for  finding 
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the  tip  diffraction.  This  will  be  done  by  using  the  6TD-MM  formulation 
and  the  magnetic  integral  equation. 

For  bistatic  scattering  a modification  to  the  method  in  this 
report  for  the  scattered  field  could  be  obtained  by  adding  the  tip  dif- 
fracted field  as  it  is  given  by  Keller,  et  al  [4].  This  formula  gives 
the  diffraction  coefficient  of  a plane  corner  with  two  straight  edges 
meeting  at  an  angle  and  is  of  the  form 


Q = _J (coso-t-cosQ*  )sin6 /24) 

4)ikQ  (cosa-cosa*  )(cosb-cosb')  ' ' 

where  a and  b are  the  angles  between  the  incident  field  and  the  two 
edges  at  the  corner,  a',  b'  are  the  angles  between  the  diffracted  field 
and  the  two  edges;  o and  o'  are  the  angles  between  the  normal  to  the 
plane  of  the  corner  and  the  incident  and  diffracted  rays;  and  6 is  the 
angle  between  the  two  edges.  Formula  (24)  along  with  expression  (22) 
will  give  at  least  the  far  zone  scattered  field.  It  will  be  seen  in 
the  next  section  that  inclusion  of  the  tip  diffracted  field  does  notice 
ably  raise  the  level  of  the  radar  cross  section  results  at  most  aspect 
angles  but  the  "scattering  pattern"  is  little  changed  by  the  tip  dif- 
fracted field. 


IV.  EXAMPLES 

In  all  examples  which  follow  we  used  a pyramid  with  the  edges  1, 

2,3  of  length  9.144X  making  an  angle  with  the  z-axis  of  15’.  At  first 
the  surface  current  was  computed.  In  this  current  there  is  no  diffracted 
current  from  the  tips.  The  scattered  field  was  computed  two  times,  one 
without  the  tip  diffraction  and  one  with.  By  using  Keller's  formula  the 
results  as  we  will  see  are  in  closer  agreement  with  the  measurements 
than  without  the  tip  diffraction. 

In  Figures  8 and  9 we  show  the  current  on  the  faces  of  the  pyramid 
(except  the  base)  for  an  incident  plane  wave  normal  to  edge  1 and  with 
the  magnetic  field  parallel  to  this  edge.  The  current  which  is  plotted 
shows  only  the  diffracted  current  since  the  current  of  the  geometrical 
optics  field  has  been  deleted  from  the  curve. 

In  Figures  10,  11  and  12  we  plotted  the  magnitude  of  the  surface 
current  of  the  pyramid  for  an  incident  plane  wave  normal  to  edge  4 and 
with  the  magnetic  field  parallel  to  this  edge.  The  bistatic  normalized 
electric  field  on  the  xy  plane  can  be  seen  in  Figure  13  for  this  case. 

In  Figures  14,  15,  16  and  17  we  show  the  bistatic  radar  cross 
section.  In  all  cases  we  keep  the  bi static  angle  constant  and  rotate 
the  pyramid  about  one  axis.  Figure  14  shows  the  bistatic  cross  section 
for  the  case  of  rotation  about  the  z axis,  with  the  transmitted  signal 
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Figure  11.  Magnitude  of 


FACE  (456) 


Figure  12.  Magnitude  of  the  surface  current 
on  face  456. 


horizontally  polarized.  In  Figure  15  we  see  the  same  situation  but  with 
the  signal  vertically  polarized.  The  bistatic  cross  sections  for  rota- 
tions about  the  y-  and  x-axes  can  respectively  be  seen  in  Figures  16 
and  17.  In  both  cases  the  polarization  is  vertical. 

As  we  can  see  in  all  cases,  by  including  the  tip  diffraction  we 
have  results  that  are  more  close  to  the  measurements  although  the 
theoretical  results  are  a little  down  from  the  measurements  probably 
because  second  order  diffraction  is  ignored.  For  the  problem  of  finding 
the  near  zone  field,  our  analysis  shows  that  we  must  take  Into  account 
diffracted  rays  of  more  than  the  first  order. 
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Figure  15.  Bistatic  (30°)  radar  cross  section  for  z polarization. 
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rotation  of  the  pyramid  about  the  x axis. 
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SUMMARY  AND  CONCLUSIONS 


The  work  reported  in  this  report  has  used  the  method  of  moments 
to  extend  the  geometrical  theory  of  diffraction  such  that  a good  approxi- 
mate solution  to  the  problem  of  scattering  by  a pyramidal  cone  is  ob- 
tained. The  solution  basically  neglects  tip  diffraction  but  it  is  shown 
that  this  effect  can  be  included  post  facto  and  improves  the  results  at 
most  aspect  angles  for  bistatic  (15°)  radar  scattering.  Even  though  the 
solution  neglects  second  and  higher  order  diffraction,  the  radar  scat- 
tering results  in  all  cases  considered  are  quite  good  as  evidenced  by 
the  agreement  between  theoretical  and  experimental  results. 

A primary  motivation  for  this  work  was  to  apply  the  moment  method 
- GTD  technique  to  a three-dimensional  geometry  to  determine  the  diffi- 
culties in  doing  so  and  also  to  determine  the  limitations  of  the  method 
when  applied  to  a three-dimensional  problem.  Thus,  while  the  formula- 
tion in  this  report  for  the  three-dimensional  problem  is  obviously  more 
involved  than  that  for  two-dimensional  problems,  it  is  not  unduly  so. 

And,  while  the  radar  cross  section  is  not  the  most  sensitive  indicator 
of  solution  accuracy,  the  results  for  the  radar  cross  section  do  indicate 
that  all  the  features  of  the  pyramid  as  a scatterer  are  well -predicted 
by  the  theory. 

As  a result  of  the  work  presented  here,  the  application  of  the 
geometrical  theory  as  extended  by  the  method  of  moments  should  find  use 
by  other  researchers  seeking  solutions  to  other  three-dimensional  problems. 
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